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Abstract:  The  paper  J'Bfuit  d’un  profil  dans  un  e<x)ulement",  by  P.  Garcia  and  P.  Gerard, 

has  been  translated  into  English.  The  paper  originally  was  published  in  Revue 
d’Acoustique,  No.  73,  1985  2/4  by  Office  National  d’Etudes  et  de  Recherches 
Aerospatiales  (ONERA).  The  work,  "Airfoil  noise  in  air  flow",  describes 
measurements  of  the  wave  number  spectrum  of  the  turbulent  boundary  layer 
induced  pressure  field  in  the  vicinity  of  the  airfoil  trailing  edge.  An  estimate  of 
radiated  noise  using  modifications  to  the  theories  of  Howe  and  Hogdson  is  made. 
The  paper  is  significant  because  of  the  good  agreement  between  the 
measurements  of  Brooks  and  Hogdson  (1980)  and  a  predicted  radiated  noise  level 
based  on  a  measured  trailing  edge  fluctuating  pressure  wave  number  spectrum 
with  measured  far  field  radiated  noise.  An  effect  of  array  element  separation 
distance  on  measured  convection  velocity  is  also  shown. 
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ABSTRACT 

An  experimental  investigation  of  the  NACA  0012  airfoil  noise  in  a  uniform  flow  has  been  conducted  in 
the  anechoic  open  flow  facility  CEPRA  19. 

The  purpose  of  these  studies  is  an  estimation  of  the  radiated  noise  deduced  from  theoretical  models 
developed  by  Howe  and  Hogdson.  These  models  need  the  spectral  characteristics  of  the  turbulent 
boundary  layer  induced  pressure  field  in  the  vicinity  of  the  trailing  edge. 

The  airfoil  has  been  instrumented  by  longitudinal  and  transverse  arrays  of  surface  pressure  sensors 
allowing  the  determination  of  the  wavenumber  spectrum  of  the  turbulent  boundary  layer.  These 
measurements  are  used  as  the  input  data  for  a  prediction  method. 
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Longitudinal  and  transverse  turbulence  integral  scales 
Airfoil  span  (1.60  m)  interacting  with  a  uniform  flow 
Airfoil  span  of  Brooks  and  Hogdson  [11]  (40  cm) 

Time  dependent  pressure  at  point  (X1.X3) 

Space-time  Fourier  transform  of  p(x,,X3,t) 

Distance  between  the  source  and  the  receiver  (1.80  m) 

Distance  between  the  source  and  the  receiver  for  Brooks  and  Hogdson  (1.20  m) 

Reynolds  number 

Effective  transducer  radius 

Acoustic  power  spectral  density 

Time 

Freestream  velocity 
Mean  convection  velocity 

Longitudinal  and  transverse  coordinates  on  the  surface  of  the  airfoil 
Coherence  function 

Displacement  thickness  of  the  turbulent  boundary  layer 
Longitudinal  and  transverse  exponential  decay  factors 

Observation  angle  in  the  horizontal  plane  (0  =  0  corresponds  to  the  downstream 
direction) 

Wave  number  spectrum  of  the  induced  fluctuating  pressure  under  the  turbulent 
boundary  layer 

Wave  number  spectrum  of  the  induced  fluctuating  pressure  in  the  proximity  of 
the  trailing  edge 

Longitudinal  and  transverse  wave  number  spectra 
Fluctuating  pressure  spectral  density 

Longitudinal  and  transverse  separations  between  two  transducers 
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Ax,,  Axj  Longitudinal  and  transverse  separations  between  pressure  transducers  that  form 

the  array 

Longitudinal  array,  short  Ax,  =  6  mm,  medium  Ax,  =  18  mm 
Transverse  array,  short  Ax,  =  3  mm,  long  Axj  =  12  mm 


p  Density 

Top,  Delay  due  to  convection,  corresponding  to  the  maximum  of  the  cross-correlation 


1.  INTRODUCTION 

The  preliminary  results  presented  in  this  article  are  relative  to  the  trailing  edge  noise  of  a  NACA  0012 
airfoil  in  a  uniform  flow,  at  low  Mach  number  (M=0.2)  and  at  Reynolds  number.  Re  =  2  x  10*, 
calculated  based  upon  a  chord  of  50  cm.  The  tests  took  place  in  the  open  jet  anechoic  wind  tunnel, 
CEPRA  19. 

The  models  developed  by  Ffowcs  Williams  and  Hall  [1],  among  others,  through  the  acoustic  analogy  of 
Lighthill;  by  Chase  [2]  and  Chandiramani  [3]  for  hydroacoustic  techniques  and,  more  recently  by  Howe 
[4]  aided  by  a  unified  theory,  attempt  to  estimate  Ae  noise  radiated  by  a  half-plane  in  a  uniform  flow. 
Howe  has  shown,  in  particular,  that  for  low  Mach  numbers  these  theories  lead  to  similar  results. 

In  this  article,  the  estimation  of  radiated  noise  by  a  NACA  0012  airfoil  will  be  made  using  a  simplified 
formulation,  given  by  Howe  [4].  Its  usage  requires  knowledge,  in  the  proximity  of  the  trailing  edge,  of 
the  wall  pressure  spectral  densities  induced  by  a  turbulent  boundary  layer  as  well  as  the  transverse 
integral  scales  of  turbulence,  iy 

These  characteristics  of  the  turbulence  are  obtained  from  experimental  results  for  a  fully  developed 
turbulent  boundary  layer  in  the  following  conditions:  an  airfoil  with  tripped  transition,  an  incidence  angle 
of  0.4' ,  and  a  velocity  of  60  m/s. 

The  experimental  model  and  instrumentation  set-up  placed  in  CEPRA  19  is  presented  in  the  first  part. 
Then  experimental  results  relative  to  the  fluctuating  pressure  field  are  given  in  the  form  of  spectral 
densities,  according  to  the  model  of  Corcos  [5],  and  wave  number  spectra. 

Finally,  the  noise  radiated  by  an  airfoil  is  identified  and  localized  due  to  the  technique  of  acoustical 
imagery  and  measured  cross-spectra,  then  compared  to  the  results  of  the  simplified  prediction 
calculations. 
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2.  TEST  CONDITIONS  AND  SIGNAL  PROCESSING  TECHNIQUES 

2.1.  Open  Flow  Facility  CEPRA  19 

The  wind  tunnel  CEPRA  19  is  an  anechoic  open  jet  wind  tunnel.  It  has  been  configured  with  a 
removable  nozzle  2.0  m  in  diameter.  Under  these  conditions,  the  test  Jet  occupies  an  approximately 
cylindrical  volume  of  which  the  diameter  and  the  length  are  2.0  and  10.6  m  respectively.  The  velocity 
of  the  flow  is  kept  between  40  and  100  m/sec. 

Figure  1  represents,  schematically,  the  installation  of  the  experimental  device  placed  in  the  anechoi?  wind 
tunnel  in  order  to  conduct  the  tests.  It  includes,  essentially,  the  NACA  0012  airfoil  and  an  array  ui 
microphones  situated  vertically  external  to  the  flow  (Figure  2). 

2.2.  NACA  0012  Model 

2.2.1.  Assembly 

The  airfoil  is  placed  vertically  in  the  wind  tunnel.  The  trailing  edge  is  situated  about  2  m  from  the  exit 
plane  of  the  nozzle  (Figures  I  and  2).  A  device  permits  adjustment  of  the  angie  of  attack  of  the  airfoil. 

The  results  presented  in  this  article  have  been  obtained  with  an  incidence  angle  near  O',  more  exactly, 
0.4',  with  the  trailing  edge  slightly  directed  toward  the  side  wall  of  the  wind  tunnel  (Figure  1). 

The  model,  of  which  the  span  and  chord  are  equal  at  3.5  ra,  is  made  up  of  a  stratified  skin,  reinforced 
by  a  metallic  frame.  A  foam  is  injected  between  the  frame  and  the  skin  in  order  to  deaden  vibrations. 
Strips  of  roughness,  2  cm  wide,  have  been  placed  at  15%  of  the  chord  from  the  leading  edge. 

2.2.2.  Instrumentation 

It  is  in  the  form  of  an  array  of  fluctuating  pressure  transducers,  flush  mounted,  and  of  static  pressure 
ports  distributed  along  the  chord  on  the  lower  face  and  the  upper  face. 

2.2.2. 1.  Arrays  (Figure  3) 

The  model  is  equipped  with  unsteady  pressure  sensors.  These  sensors  are  distributed  in  arrays:  medium, 
short  and  long,  and  longitudinal  and  transverse  arrays.  This  instrumentation  is  situated  near  the  trailing 
edge.  Two  types  of  sensors  are  used: 

*  Entran  EPIL  203A  -  They  have  a  diameter  of  2.03  mm  and  of  a  length  of  5  mm,  and  are 
integrated  in  a  small  transverse  bar  of  about  100  mm  within  3  mm  of  the  trailing  edge.  There  are  15  of 
them.  These  sensors  are  placed  near  the  trailing  edge  parallel  to  the  surface  of  the  airfoil,  and  connected 
to  the  latter  by  a  cylindrical  cavity  of  approximate  length  2  mm  and  diameter  1  mm.  A  study  of  the 
amplitude  and  phase  response  of  the  transducer  in  the  presence  of  this  cavity  has  been  made  by 
CERT/ONERA  [6].  The  acoustical  resonance  frequency  has  been  identified  as  13  kHz.  Accounting  for 
the  diameter  of  the  cavity,  the  cutoff  frequency  of  the  device  is  estimated  [7],  for  an  attenuation  of  5  dB 
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Figure  1.  Experimental  set-up  of  the  model  and  the  vertical  antenna  in  CEPRA  19. 

^  micro  79 


Figure  2.  Experimental  set-up  of  the  NACA  0012  model  and  the  instrumentation. 
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Direction  of  Flow 


Figure  3.  Trailing  edge  instrumentation, 
of  the  wall  pressure,  approximately  by: 

0.8  X  60 

f  =  -  =  15  kHz 

2ir  X  0.5  X  10'^ 

*  Electret  sensors  -  These  sensors  are  based  on  the  "electret"  effect,  which  is  known  as  the  quasi¬ 
permanent  electrostatic  polarization  of  a  plastic  membrane,  of  which  the  geometric  radius  is  1.8  mm. 
In  these  conditions  of  measuring  wall  pressure,  its  effective  radius  is  about  1.08  mm  (0.6  x  1.8)  where 
the  cutoff  frequency  is  estimated  to  be  7  kHz.  They  are  placed  in  line  along  the  chord,  perpendicular 
to  the  bar  of  Entran  sensors. 

2. 2. 2. 2.  Static  Pressure  Ports 

A  series  of  47  local  static  pressure  taps  have  been  arranged  along  the  chord  on  the  lower  and  upper  faces 
in  order  to  measure  the  pressure  distribution.  These  taps  are  linked  by  a  pressure  switch  to  48  supply 
lines  from  a  differential  pressure  transducer  rated  at  ±34.5  kPa.  The  pressure  distribution  along  the 
chord  allows  the  calculation  of  the  incidence  of  the  airfoil. 
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2.3.  Instrumentation  outside  of  the  model 

It  is  made  up  of  a  multimicrophone  array  (acoustic  antenna)  and  of  2  B&K  model  4135  microphones 
placed  symmetrically.  The  multimicrophone  array  is  composed  of  28  B&K  microphones  spaced  10  cm 
apart,  on  a  2.7  m  long  beam  placed  outside  of  the  flow  and  parallel  to  the  trailing  edge  of  the  airfoil 
(Figure  2). 

2.4.  Acquisition  and  signal  processing 

For  each  test,  the  signals  from  the  different  transducers  are  amplified  (research  amplifiers  that 
automatically  set  the  gain,  developed  by  ONERA)  and  are  recorded  on  three  28  track  tape  recording 
devices.  Parallel  to  these  recording  devices,  a  Fourier  analyzer,  model  HP  5451  B  dealt  with,  in  real 
time,  the  signals  of  any  three  transducers  allowing  initial  results  to  be  acquired  during  the  test.  The  set 
of  results  obtained  in  the  Fourier  system  is  stored  numerically  on  magnetic  tape  in  order  to  be  played 
back  on  the  CYBER  system  at  ONERA.  A  complete  package  for  analyzing  the  measurements  and  for 
graphical  output  has  been  developed  at  ONERA. 

The  multisensor  computations  are  made  with  the  aid  of  a  vector  processor  AP  120  (Floating  Point 
System),  connected  to  an  HP  1000  computer.  In  this  way,  one  obtains  an  acoustic  image  of  sources  of 
noise.  The  adaptation  of  this  technique  in  the  near  field  for  tests  in  the  wind  tunnel  takes  into  account 
the  effects  of  convection  and  refraction  and  is  discussed  in  [8]. 

Two  types  of  analysis  exist.  The  first  is  analog  and  is  obtained  in  real  time,  from  which  it  get  its  name, 
"Real  Time  Antenna".'  It  consists  of  summing  the  delayed  signals  of  each  of  the  microphones  of  the 
antenna,  for  a  given  focal  point,  and  establishing  the  spectrum  of  the  summed  signal.  The  second  is 
numerical  and  i:  cbrair.;  J  off  line,  it  givvs  on  acoustic  Image  that  plots  KX)  spectra  corresponding  to  100 
focal  points  (Figure  14). 

Similarly,  the  programs  to  calculate  the  wave  number  spectra  have  been  developed  at  ONERA.  A 
completely  similar  analysis  is  given  in  reference  [9]. 


'  [Translators’  note:  Ther  term  "antenna"  is  meant  in  the  sense  of  capable  of  being  directive.  In  the 
United  States,  this  would  be  called  a  beamformable  or  directional  array.] 
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3.  FLUCTUATING  PRESSURE  HELD  INDUCED  BY  A  FUTLY  DEVELOPED  TURBLXENT 
BOUNDARY  LAYER 

Preliminary  measurements  have  been  made  to  characterize  the  fluctuating  pressure  field  induced  by  a  fully 
developed  turbulent  boundary  layer  at  the  velocity  of  60  ni/'sec. 

Tne  results  presented  were  obtained  from  signals  furnished  by  the  transducers  comprising  the  longitudinal 
and  transverse  arrays. 

The  two  types  of  analysis  which  have  been  used  to  characterize  the  fluctuating  wall  pressures  are  based 
one.  on  the  longitudinal  and  transverse  cross-spectral  density  functions  (Corcos  model  [5]).  and  the  other, 
on  the  wave  number  spectra. 

3.1.  Theoretical  model,  Corcos  [5],  and  turbulence  integral  scales 

The  arrays  of  pressure  transducers  are  found  near  the  trailing  edge,  in  a  region  of  a  weak  positive 
pressure  gradient.  Under  these  conditions,  the  Corcos  model  [5],  developed  for  a  flat  plate,  remains 
valid.  This  has  been  verified  experimentally  by  Schloemer  (10). 

The  Corcos  model  describes  a  pressure  field  induced  on  a  surface  by  a  homogeneous  turbulent  boundary 
layer.  The  two  dimensional  field  is  expressed  separately  according  to  the  longitudinal  and  transverse 
components: 

(1)  'J’TBL(ei.€3.«)  =  A(wei/V)  BCwCj/V)  exp(-ia)£i/V) 

where  e,,  e,  are  the  longitudinal  and  transverse  separations  for  two  points  on  the  surface  of  ^e  airfoil. 
V  is  the  convection  velocity,  is  the  pressure  spectral  density  and  A(«€,/V)  and  Blwcj/V)  are  the 
longitudinal  and  transverse  Corcos  functions. 

For  the  longitudinal  and  transverse  directions,  the  turbulence  integral  scales  are  defined  in  the  following 
way  [11]: 

'■  ■  (  #)  '■  '  \  "*' 

and 

With  the_fonctions  A(«Ci/V)  =  exp(-7;,  «e,/V)  and  B(wcj/V)  =  exp(-T;j  ue^lV)  and  a  constant  convection 
velocity  V,  the  previous  equations  are  simplified  [3,11]  to  give  the  following  expressions: 


(4) 


where  tj,  and  are  the  decay  factors  of  the  functions  A((i;€,/V)  and  BioiCylV). 
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The  Corcos  functions  A(wei/V)  and  B(a)€,/V)  correspond  to  cross-spectral  functions  of  which  the  abscissa 
is  put  into  a  nondimensionaJ  form.  They  are  derived  from  the  classic  coherence  of  the  type 
7  =  /  V  .  where  G^  and  G„  are  the  autospectra  of  the  signals  u  and  v  and  G„^  is  the  cross¬ 

spectrum  of  these  signals. 

This  definition,  slightly  different  from  that  used  by  Corcos,  nevertheless  remains  valid  where  the 
assumption  of  homogeneity  is  respected. 

The  coherences  are  calculated  from  signals  furnished  by  the  transducers  located  in  provimity  to  the 
trailing  edge. 

For  the  longitudinal  array,  transducer  number  16  has  been  chosen  as  the  reference  (Figure  3).  Under 
these  conditions,  the  separation  distances  between  the  transducers  vary  from  6  to  72  mm.  Figure  4 
represents  the  development  of  the  longitudinal  coherences,  in  the  frequency  band  between  I  and  10  kHz. 

The  longitudinal  coherences  give  rise  to  two  remarks: 

-  the  curves  show  a  plateau  at  about  1  kHz,  no  matter  what  the  gap  between  the  transducers; 

-  beyond  1  kHz,  the  curves  decrease  rapidly. 

For  the  transverse  array,  transducer  number  7  has  been  chosen  as  the  reference  (Pigure  3). 

Although  the  transverse  distances  are  shorter,  3  to  15  tnm,  the  coherence  functions  are  significant  over 
a  much  more  narrow  frequency  band  compared  to  the  longitudinal  ones  (Figure  5). 

These  transverse  coherences  are  marked  by  a  distinctly  pronounced  hump  whose  maximum  is  located  near 
1  kHz,  no  matter  what  the  gap;  beyond  this  value,  these  curves  decrease  rapidly. 


Figure  4.  Longitudinal  coherence  functions  analyzed  up  to  10  kHz,  for  different  separations  e,  - 
Transducer  16  as  reference  -  U,  =  60  m/sec. 
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Figure  5.  Transverse  coherence  functions  analyzed  up  to  10  kHz,  for  different  separations  e,  - 
Transducer  7  as  reference  -  =  60  m/sec. 

To  transform  the  abscissas  of  these  nondimensional  curves,  in  the  form  we/V,  requires  the  choice  of  a 
convection  velocity  V.  Figure  6  groups  together  the  normalized  cross-correlations  between  the  signals 
coming  from  transducer  16,  the  reference,  and  the  sensors  which  comprise  the  short  and  medium 
longitudinal  arrays.  The  propagation  times  are  then  presented  on  Figure  7  as  a  function  of  longimdinal 
separations  The  convection  velocity  is  calculated  from  the  slope  of  the  curve  (Tigure  7)  and  is  close 
to  0.7  Uo,  when  the  separations  referenced  to  transducer  16,  are  less  than  40  mm.  A  higher  velocity. 
0.8  U„,  is  obtained  for  larger  separations  e,.  This  latter  one  has  been  kept  to  nondimensional ize  the 
longitudinal  and  transverse  coherence  functions. 

Figure  8  and  9  plot  these  results  inferred  from  Figures  4 _Md  5  for  different  separations  e,  and  ci.  From 
these  sets  of  points  it  is  clear  the  average  curves  A(«e,/V)  and  BfucjAO  conform  to  the  Corcos  model. 
The  function  Bfucj/V)  decreases  much  more  rapidly  than  A(we,/V).  These  functions  are  described 
analytically  by  decreasing  exponentials  of  which  the  decay  factors  are  equal  to  0.17  and  0.9  respectively 
according  to  the  longitudinal  and  transverse  coordinates. 
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I  Normalized  Cross-correlation 


I 


Figure  6.  Normalized  cross-correlations  for  different  length  longitudinal  arrays  of  transducers  - 
Uq  =  60  m/sec,  —  short  array  Axj  =  6  mm, —  medium  Axj  =  18  mm 


Figure  7.  Variation  of  rgpj  as  a  function  of  longitudinal  separations  -  Uq  =  60  m/sec, 
—  short  array  AX|  =  6  mm,  — medium  Axj  =  18  mm 
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Figure  8.  Longitudinal  coherence  as  a  function  of  a  nondimensional  number  u£j/V.  Uq  =  60  m/sec. 


Figure  9.  Transverse  coherence  as  a  function  of  a  nondimensional  number  we^/V.  Uq  =  60  m/sec. 

For  a  frequency  of  1000  Hz,  the  turbulence  integral  scales,  given  by  equation  (4),  give  the  following 
approximate  values: 

0.8x60 

= -  =  45  X  la’  ra  (1.77  in.) 

2tx0.17x  10* 


Airfoil  Noise  in  an  Air  Flow 


13 


with  T/i  =  0. 17 

0.8  X  60 

and  =  -  =  8.5  x  10"^  m  (0.33  in.) 

It  X  0.9  X  10^ 

with  Tjj  =  0.9. 

Table  1  shows  similar  results,  obtained  by  Hogdson  and  Brooks  [11],  for  a  NACA  0012  airfoil  with  a 
chord  of  60  cm.  It  is  noted  that  the  transverse  decay  factor  jjy  of  0.9  as  well  as  the  convection  velocity 
of  0.8  Uo  are  greater  than  those  obtained  by  Hogdson. 

3.3.  Theory  and  definitions;  wave  number  spectrum 

The  fluctuating  pressure  field  induced  by  a  fully  turbulent  boundary  layer,  has  been  described,  until  now, 
in  terms  of  spatial  cross-spectral  functions  according  to  the  theory  of  Corcos. 

Another  way  of  modeling  the  homogeneous  hydrodynamic  pressure  field  is  to  consider  it  as  a  distribution 
of  evanescent  harmonic  waves  (3).  The  pressure  wave  number  spectrum  is  derived  from  the  double 
space-time  Fourier  transform  of  the  wall  pressure  p(x„X3,t),  at  the  point  Xj.x,  taken  on  the  surface  of  the 
airfoil. 

F.T. 

(5)  p(x,,X3,t) - »  P(ki,k3,«) - ^  $xBL(k„k3,w)  =  1  P(k„k3,u)  |  ^ 

In  using  the  Corcos  hypothesis,  the  expression  for  the  wave  number  spectrum  of  the  hydrodynamic 
pressure  field  simplifies: 

(6)  #TBL(k„k3,u)  =  'l'(aj)  f-Ik,)  ^>(k3) 

where  the  functions  $(w),  ‘J'Ck,)  and  ^fk,)  are,  respectively,  the  fluctuating  pressure  spectral  density  and 
the  spatial  Fourier  transforms  of  the  Corcos  functions. 


lA  V 


-iwe,  t 

exp  exp(-ik,e,)  d«, 


with 
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and 

(8) 


<^3  N 


2n 


\ 


} 


exp(-ikj  fij)  de, 


with 


for  a  constant  convection  velocity  V,  the  spectral  functions  $i(ki)  and  'J’jCkj)  are  simplified  [3]: 


(9) 


and 


(10) 


3.4.  Experimental  results  for  wave  number  spectra 

The  wave  number  spectra  •J’lClTi.f)  and  ^jCkj.O,  presented  in  Figures  10  and  12,  are  obtained  from  the 
short  longitudinal  and  short  transverse  arrays,  for  which  equidistant  distributions  of  Ax,  and  Ax,  are  6 
and  3  mm,  respectively. 

These  results  are  presented  in  a  three  dimensional  form  with  the  amplitudes  4'i(ki,f)  and  ^sCki.f), 
expressed  as  a  function  of  k,,  Icj,  and  f,  where  ki  and  ks  are  defined  as  being  the  ratio  of  a  frequency  over 
a  velocity. 

Figure  10  shows  the  appearance  of  spectral  aliasing  in  k  for  frequencies  above  2500  Hz.  In  this  example, 
obtained  with  a  spatid  sample  Ax,  of  6  mm,  the  aliasing  in  k  should  be  produced  whenever  k  ^ 

=  1/2  Ax,  =  83  m  ‘.  This  value  is  effectively  reached  by  the  components  of  the  hydrodynamic  pressure 
field,  convected  at  the  velocity  of  0.7  U,  (Figure  11),  when  the  frequencies  are  equal  to  or  greater  than 
2500  Hz. 

The  energy  of  $,(r,,0  is  concentrated  in  the  convective  region  in  which  the  trace  in  the  plane  (ki,f),  at 
the  mean  of  the  grey  shaded  representation,  gives  approximately  a  straight  line  (Figure  11). 

The  convective  velocity,  derived  from  the  slope  of  this  straight  line,  is  equal  to  0.7  U,.  A  similar  result 
for  $,(ir,,f)  has  been  investigated  from  the  medium  longitudinal  array,  where  the  separation  Ax,  between 
the  transducers  is  equal  to  18  mm.  The  preceding  remarks  remain  true,  but  the  velocities  obtained  are 
equal  to  0.8  U,.  These  observations  parallel  the  conclusions  drawn  from  Figure  7. 

It  follows  then  that  the  medium  and  small  scale  turbulent  boundary  layer  structures  are  spatially  filtered 
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by  the  longer  arrays  and  that  their  convection  velocity  0.7  U,  is  less  than  those  of  structures  of  much 
larger  dimensions  which  stabilize,  from  the  array  being  longer,  in  the  vicinity  of  0.8  U„. 

Distance  between  sensors:  6.0  mm 
500  Averages 
Point  Number:  22 
Number  of  sensors:  10 
U,:  60.0 


0,  (R 


0.00252 
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0.00105 
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0 


Wave  Number  Spectrum  (Kj,F) 


Figure  10.  Wave  number  spectrum  <J>i(k„f)  obtained  from  the  short  longitudinal  array.  U„  =  60  m/sec. 
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Figure  11.  Projection  of  the  wave  number  spectrum  on  the  (k„f)  plane. 
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These  experimental  values  of  convection  velocity  remain  greater  than  those  of  Hodgson  (see  Table  1). 


Table  1. 


Study 

’ll  (m) 

’?3  (m) 

V 

CONDITIONS 

Garcia  -  Gerard 

+  0.17 

+  0.9 

-  0.8  U. 

12%  NACA  Profile 
Chord  50  cm 

U„  =  uO  m/sec 

Hogdson  -  Brooks 

+  0.14 

+0.58 

=  0.8  U„ 

12%  NACA  Profile 
Chord  60  cm 

Uo  =  69  m/sec 

Corcos 

+0.11 

+0.71 

Flat  plate 

U„  =  69  m/sec 

The  wave  number  spectrum  ^jCkj.f),  obtained  from  the  transverse  array,  is  symmetrical  in  relation  to  the 
plane  k,  =  0  and  centered  around  the  frequency  of  I  kHz  (Figure  12). 
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Figure  12.  Wave  number  spectrum  'f3(E3,0  obtained  from  the  short  transverse  array.  U,  =  60  m/sec. 
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4.  RADIATED  NOISE 

The  prediction  calculation,  based  on  the  theory  of  Howe  [4],  is  given  in  the  beginning  of  this  paragraph. 

The  noise  radiated  by  the  airfoil  is  localized  and  identified  by  acoustic  imagery  techniques  and  from 
analysis  of  cross  spectra.  The  measurement  of  this  noise,  due  to  the  interaction  of  the  stable  flow  with 
the  airfoil,  is  obtained  beginning  with  the  autospectra  of  signals  furnished  by  two  microphones,  placed 
symmetrically.  A  comparison  between  the  measured  and  calculated  radiated  noise  from  a  simplified 
formulation,  concludes  this  chapter. 

4.1.  Calculated  predictions  based  on  the  theory  of  Howe 

Howe  [4]  proposes  an  extension  of  Lighthill’s  analogy  to  a  stationary,  irrotational,  non-viscous,  non¬ 
conducting,  and  isentropic  mean  flow.  The  solving  of  the  wave  equation  for  enthalpy  allows  one  to 
express: 


-  the  "hydrodynamic"  pressure  induced  by  the  turbulent  boundary  layer,  in  the  neighborhood 

of  the  trailing  edge;  and 

-  the  radiated  acoustic  pressure  in  the  far  field. 

Howe  obtained  a  relation  between  the  radiated  noise  spectrum  and  the  wave  number  spectrum 
1'TE(k,,  w  cosa/V,  u)  of  the  hydrodynamic  pressure  that,  in  the  far  field,  is  written; 


where  V  equaling  0.8  U„,  is  the  convection  velocity,  L,  the  span,  R,  the  distance  from  the  trailing  edge 
to  the  receiver,  and  6,  the  angle  of  observation. 

For  an  observation  angle  of  ^  =  90' ,  the  preceding  equation  reduces,  according  to  Brooks  and  Hodgson 
til],  to. 


(12) 


S(«) 


1 

>/2  X*  R* 


♦(«) 


with  y<I'TE(k„0,aj)  dk,  =  fjCw)  $(u). 


4.2.  Radiated  Noise  Measurement 

4.2.1.  Analysis  of  radiated  noise  using  microphone  63  and  acoustic  imagery. 

The  first  estimate  of  radiated  noise  by  the  airfoil  is  obtained  by  direct  comparison  of  the  spectra  of  the 
acoustic  pressure  emitted  by  the  airfoil  to  those  of  the  background  noise  of  the  wind  tunnel. 
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Figure  13  presents  tiie  spectral  analyses  made  in  the  frequency  band  300  Hz  -  5  kHz,  from  signals 
furnished  by  microphone  63.  The  following  points  are  emphasized: 

-  the  emission  frequency  band  for  the  airfoil  is  located,  in  this  first  approach,  in  a  few  hundred 

hertz  around  2  kHz;  and 

-  comparative  study  of  the  spectra  (Figure  13)  makes  apparent  two  very  distinct  regions:  one 

corresponds  to  spreads  greater  than  10  dB,  high  pass  filtered  at  300  Hz,  in  the  frequency 
band  less  than  1  kHz,  the  other  to  spreads  greater  than  5  dB  for  medium  and  high 
frequencies. 

Figure  14  represents  an  acoustic  image  composed  of  100  focal  points  distributed  symmetrically  along  the 
length  of  the  trailing  edge  over  a  distance  of  +  2.50  m.  When  the  airfoil  is  present  in  the  test  flow,  two 
sources  stand  out  distinctly.  The  acoustic  imagery  sections  (Figure  15),  for  frequencies  from  850  to  1000 
Hz,  localize  a  source  at  1.15  m  and  the  other  at-0.68  m. 

In  effect,  the  airfoil  with  a  3  m  span  interacts  with  a  jet,  shown  schematically  on  Figure  2,  composed  of 
a  potential  cone  and  a  mixing  zone.  These  two  types  of  flow,  one  laminar  and  one  turbulent,  give  birth 
to  two  sources,  a  low  frequency  one  in  the  mixing  zone,  the  other  one  at  medium  and  high  frequencies 
in  the  potential  cone. 

The  two  parasitic  sources  localized  by  the  acoustic  imagery  technique  (Figure  14)  result  from  the 
interaction  of  the  mixing  zone  with  the  airfoil. 

4.2.2.  Localization  and  identification  of  the  acoustic  source,  coming  from  interaction 
with  the  potential  cone 

The  localization  and  the  identification  of  the  sought  after  aerodynamic  source  has  been  established  from 
cross-spectral  analysis  between  pressure  transducers  near  the  trailing  edge  and  a  "far  field"  receiver. 

In  order  to  partly  separate  itself  from  these  two  interfering  sources,  the  "real  time"  acoustic  antenna  has 
been  focusal,  using  analog  methods,  on  the  instrumented  part  of  the  airfoil. 

The  normalized  cross-correlation  function  (Figure  16)  between  the  analog  signals  furnished  by  the  "real 
time"  antenna  and  pressure  transducer  7,  analyzed  from  500  Hz  to  2  kHz,  shows  a  delay  of  14.5  ms, 
which  effectively  corresponds  to  the  time  taken  by  the  sound  pressure  to  propagate  from  a  region  situated 
near  the  trailing  edge  towards  the  antenna. 

We  emphasize  that,  for  the  relative  positions  of  the  source  and  the  receiver  in  relation  to  the  mixing  zone 
of  the  jet,  the  combined  influence  of  the  effects  of  convection  and  refraction  on  the  propagation  is 
negligible  [12]. 

The  corresponding  cross-spectrum  (Figure  17)  analyzed  in  a  frequency  band  of  500  Hz  -  4  kHz,  gives 
a  phase  plot  which,  although  noisy,  settles  at  a  value  of  -t/4  ^  in  the  500  Hz  to  2  kHz  range.  This  value 
of  the  phase  is  obtained  after  having  taken  into  account  the  correction  for  the  electronic  phase  differences 

^  [Translators’  note:  The  phase  plot  provided  in  Figure  17  actually  shows  this  as  +3x/4.  So  there  must 
have  been  a  -t  electronic  phase  correction.] 


Figure  13.  Comparison  of  radiated  spectra  with  and  without  airfoil  in  the  flow.  Microphone  63 
Uq  =  60  m/sec. 

.  Arbitrary  units 


Figure  14.  Soanwise  "acoustic  image"  of  the  NACA  0012  airfoil,  Uq  =  60  m/sec. 


Spanwise  location  (m) 

Figure  IS.  Sections  obtained  from  the  acoustic  image  for  frequencies  of  850  and  1000  Hz. 
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Figure  16.  Cioss-correiation  functions  between  the  acoustic  antenna  (real-time)  and  the  flush  mounted 
pressure  transducers  indicated. 

a)  transducer  7 

b)  transducer  17 

c)  transducer  21 

d)  transducer  24 
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Figure  17.  Cross-spectrum  between  the  acoustic  antenna  (real-time)  and  flush  mounted  pressure 
transducer  number  7.  Frequency  range  of  analysis  5(X)  Hz  -  4  kHz,  Uq  =  60  m/sec. 
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Figure  18.  Cross-spectrum  and  autospectrum  from  microphones  78  and  79,  Uq  -  60  m/sec. 
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of  the  Entran  transducers  with  the  B&K  microphones.  Moreover,  the  correlation  between  the  signals, 
furnished  by  transducers  17,  21,  24  (Figure  3)  and  the  antenna  becomes  even  weaker  since  the  transducer 
considered  is  displaced  from  the  trailing  edge  (Figure  16). 

These  results  permit  us  to  state  that: 

-  the  source  is  localized  in  the  region  of  the  trailing  edge;  and 

-  the  sought  after  acoustic  source  emits  in  a  frequency  band  of  about  500  Hz  -  2  kHz 

4.2.3.  Measure  of  radiated  noise  using  two  microphones  78  and  79  symmetric  about 
the  airfoil 

The  techniques  put  in  place  previously  have  allowed  the  localization  and  identification  of  radiated  noise 
by  the  trailing  edge  of  the  airfoil  situated  in  the  steady  part  of  the  flow.  On  the  other  hand,  they  do  not 
allow  us,  currently,  to  recover  the  radiated  noise  spectrum. 

For  that,  we  have  taken  the  approach  of  Brooks  and  Hogdson  [11]  that  consists  of  extracting,  from  the 
background  noise  of  the  wind  tunnel,  the  noise  emitted  by  the  airfoil  by  using  the  cross-spectrum  obtained 
from  the  signals  furnished  by  the  microphones  placed  symmetrically  apart  and  from  others  on  the  model. 

The  results,  presented  in  Figure  18,  are  obtained  for  microphones  78  and  79  (Figure  2).  The  spectral 
analysis  has  been  high  pass  filtered  above  the  frequency  of  800  Hz  in  order  to  distinguish  it  from  the 
parasitic  low  frequency. 

The  plotted  curves,  in  Figure  18,  confirm  the  selective  character,  in  the  high  frequencies,  of  the  cross¬ 
spectrum  from  the  signals  furnished  by  microphones  78  and  79,  compared  to  the  radiated  noise  spectrum. 
The  phase  of  the  cross-spectrum  is  very  well  defined  and  shows  an  opposition,  in  the  frequency  band  of 
4.5  kHz  that  characterizes  the  dipole  nature  of  the  studied  source. 

The  slope  of  these  phases  is  probably  due  to  a  very  slight  variation  between  the  propagation  times,  on 
the  order  of  0.1  ms,  of  the  acoustic  pressure  that  reaches  microphones  78-79,  for  which  the  positions  are 
not  perfectly  symmetric. 

Another  experimental  ie*jlt  of  measured  noise  is  presented  in  Figure  19  for  a  velocity  of  70  m/sec,  but 
with  the  span  and  the  distance  from  the  trailing  edge  to  the  receiver  reduced,  allowing  a  comparison  with 
the  measurements  of  Brooks  and  Hogdson  [1  Ij. 

The  reduction  factor  depends  only  on  the  geometry  of  the  airfoil,  the  spectral  densities  4'(w),  and  the 
turbulence  integral  scales  being  considered  fixed.  The  factor  is  deduced  from  the  simplified  formula 
[12]  of  Howe  and  is  equal  to  10  log  (L„/L)  +  10  log  (R/RH)^  where  Lh  and  Rh  are,  respectively,  the 
span  length  and  distance  from  the  trailing  edge  to  the  receiver  for  the  Hogdson  model  and  L  and  R 
correspond  to  those  of  this  study. 

Only  the  chord  length,  60  cm  for  Hogdson  and  50  cm  for  this  study,  differ,  A  comparison  of  this 
measurement  reduced  from  radiated  noise  with  that  of  Hogdson,  for  a  NACA  0012  airfoil  of  span  40  cm 
and  a  distance  of  1.2  m  from  the  trailing  edge  to  the  receiver  and  for  a  velocity  of  70  m/sec,  results  in 
a  overlaying  of  the  curves  (Figure  19). 
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Figure  19.  Comparison  with  the  radiated  noise  measurement  obtained  by  Brooks  and  Hogdson  (  — ), 
Uy  =  70  m/sec. 
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Figure  20.  Comparison  of  measured  and  calculated  radiated  noise.  NACA  0012  airfoil. 
Uq  =  60  m/sec. 


P.  Garcia  and  P.  Gerard 


24 


4.3.  Overlaying  of  calculated  and  measured  radiated  noise  for  a  .\ACA  0012  profile  in  a 
uniform  flow 

The  predictior.  calculation,  based  on  the  simplified  formulation  [12]  of  Howe,  is  applied  to  the 
radiated  noise  by  the  NACA  0012  airfoil,  for  the  flow  speed  Uo  =  60  m/sec.  The  calculations  require, 
in  addition  to  the  geometry  of  the  airfoil  -  receiver  group,  two  other  input  data: 

-  the  spectral  density  of  the  fluctuating  pressure  near  the  trailing  edge  (the  spectral 

density  of  transducer  17  is  represented  in  Figure  20);  and 

-  the  transverse  turbulence  integral  scale  of  fj,  deduced  from  measurements  of  the  wall 

pressure  (paragraph  3.2.). 

Only  the  interaction  between  the  potential  cone  and  the  airfoil  is  taken  into  account.  Consequently,  the 
sought  after  source  is  distributed  along  the  length  of  the  span  submitted  to  a  steady  flow,  that  is  about 
1 .6  m. 

In  Figure  20,  the  prediction  calculation  is  compared  to  the  measured  radiated  noise  by  the  airfoil, 
obtained  previously  from  the  cross-spectrum  between  the  two  microphones  78  and  79;  the  difference  is 
of  the  order  of  a  few  decibels  in  the  frequency  band  800  Hz  -  3  kHz. 

This  technique  has  been  used  by  Hogdson  [1 1  j.  The  difference  between  the  experimental  results  and  the 
calculated  results  is  analogous  to  the  one  presented  in  Figure  20. 


5.  CONCLUSIONS 

Acoustical  imaging  and  cross-spectral  techniques  allow,  on  the  one  hand,  the  highlighting  of  an  acoustic 
source  in  the  vicinity  of  the  trailing  edge  and,  on  the  other  hand,  identifying  the  range  of  frequencies  of 
emission  of  this  source  (500  Hz  -  2  kHz). 

The  two  microphones  78  and  79,  symmetrically  placed  in  relation  to  the  airfoil,  allow  the  establishment 
of  the  dipole  nature  of  the  studied  source,  as  well  as  its  measurement. 

A  comparison  of  the  measured  results  with  those  obtained  by  Brooks  and  Hogdson,  under  similar  velocity 
conditions,  similar  span,  and  trailing  edge  to  receiver  distance,  is  entirely  satisfactory. 

The  estimation  of  noise  radiated  by  the  trailing  edge  from  the  theory  of  Howe,  requires,  in  addition  to 
the  airfoil-receiver  geometry,  knowledge  of  the  wall  pressure  field,  induced  by  the  fully  developed 
turbulent  boundary  layer  established  in  the  region  of  the  trailing  edge. 

Two  techniques  are  used  in  order  to  model  this  boundary  layer: 

-  the  wave  number  spectrum  of  the  fluctuating  pressure,  ^>TE(lti.  kj,  u);  and 

-  the  pressure  spectral  density  4>(a;)  and  the  turbulence  integral  scales,  and  f,,  according  to 

the  model  of  Corcos. 
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The  input  data  #(w)  and  i,,  deduced  from  wall  pressure  measurements,  are  introduced  in  the  simplified 
expression  by  Howe  in  order  to  estimate  the  radiated  noise  in  the  far  field,  at  the  angle  6=90’ .  The 
measurement  and  the  estimation  of  radiated  noise  by  the  trailing  edge  of  an  airfoil  agrees,  within  a  few 
decibels,  in  the  frequency  range  of  800  Hz  -  3  kHz,  for  a  velocity  of  60  m/sec. 

Our  subsequent  objectives  are  to  use  the  more  general  estimation  from  Howe  who  introduced  a  model 
of  the  turbulent  boundary  layer  in  the  form  of  a  wave  number  spectrum.  The  predicted  noise  levels  will 
be  verified  by  other  flow  conditions  at  the  trailing  edge,  obtained  by  varying  the  airfoil  leading  edge. 
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